Study on the large scale dynamo transition 
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Using the magnetohydrodynamic (MHD) description, we develop a nonlinear dynamo model that 
couples the evolution of the large scale magnetic field with turbulent dynamics of the plasma at small 
scale by electromotive force (e.m.f.). The nonlinear behavior of the plasma at small scale is described 
by using a MHD shell model for velocity field and magnetic field fluctuations. The shell model allow 
to study this problem in a large parameter regime which characterizes the dynamo phenomenon in 
many natural systems and which is beyond the power of supercomputers at today. Under specific 
conditions of the plasma turbulent state, the field fluctuations at small scales are able to trigger 
the dynamo instability. We study this transition considering the stability curve which shows a 
strong decrease in the critical magnetic Reynolds number for increasing inverse magnetic Prandlt 
number Pm" 1 in the range [10~ 6 , 1] and slows an increase in the range [1, 10 s ] . We also obtain 
hysteretic behavior across the dynamo boundary reveling the subcritical nature of this transition. 
The system, undergoing this transition, can reach different dynamo regimes, depending on Reynolds 
numbers of the plasma flow. This shows the critical role that the turbulence plays in the dynamo 
phenomenon. In particular the model is able to reproduce the dynamical situation in which the 
large-scale magnetic field jumps between two states which represent the opposite polarities of the 
magnetic field, reproducing the magnetic reversals as observed in geomagnetic dynamo and in the 
VKS experiments. 
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One of the most fascinating and challenging topics in 
physics and astrophysics is the understanding of the gen- 
eration and self-sustaining of magnetic fields in planets, 
stars, galaxies, etc. The most accredited mechanism is 
the so-called dynamo effect, i.e. the maintaining of a 
magnetic field against diffusive effects by the motion of 
electrically conducting fluids. This effect has also been 
studied in laboratory liquid metal experiments like Riga 
and VKS experiment [l| and plays a fundamental role 
in our understanding of many magnetized phenomena 
which are of interest in many research fields since Fara- 
day's study. 

The dynamo effect occurs because magnetic field lines 
are generally stretched at small scales by the random mo- 
tion of the fluid in which they are almost 'frozen'. The 
magnetic lines stretching leads to an exponential amplifi- 
cation of the field, until the backreaction (via the Lorentz 
force) causes the saturation of this growth. The small- 
scale dynamo is generated when the system reaches en- 
ergy equipartition at small scales. During the growth 
of the magnetic energy at small scales, the velocity and 
magnetic field fluctuations interacts each other, generat- 
ing an e.m.f. able to produce a large scale magnetic field. 
The latter increases due to the turbulent fluctuations at 
small scales until a saturation level is reached. The dy- 
namo quenching occurs because the system attempts to 
conserve the total magnetic helicity when the magnetic 
field is growing [2|. 

In the present Letter we investigate the dynamo effect 
driven solely by turbulent fluctuations, in absence of a 
mean flow (a 2 -dynamo), focusing our study on the dy- 
namical transition towards the regime where the large 



scale magnetic field starts to be generated. Dynamo 
transition results from an instability: when, for an as- 
signed value of the Reynolds number Re ~ Su/kov, the 
magnetic Reynolds number Rm ~ Su/kofi (5u being the 
r.m.s. of the turbulent velocity, v and fi respectively the 
kinetic viscosity and the magnetic diffusivity) reaches a 
critical value Rm c , the magnetic field looses its stabil- 
ity from a quasi zero-magnetic field state generating the 
initial exponential increasing in its values until the dy- 
namo quenches. Generally, the dynamo bifurcation could 
display at least two different natures: subcritical and su- 
percritical [3j, in analogy with turbulent transitions. We 
want to investigate its nature, using a new model, re- 
cently built up, with the aim to overcome the limitation 
in the Reynolds number range covered by both direct nu- 
merical simulation (DNS) and by a lot of other theoretical 
models. In fact, although DNS are playing an important 
role in our understanding of the dynamo problem [J], a 
realistic parameter regime is still beyond the power of 
supercomputers at today. Difficulties arise from the very 
large values of the dimensionless parameters, which char- 
acterize this problem in natural systems. 

In 1955, Parker Q suggested that the net effect of av- 
eraging many small scale turbulent motions would be to 
produce the large scale electric field (a effect) generat- 
ing large scale poloidal and toroidal magnetic field. The 
latter can also be generated by the differential rotation 
(a-Q, effect). In the same spirit (see also [y]), we make 
a decomposition of the fields in an average part, vary- 
ing only on the large scale L, and a turbulent fluctuating 
part, varying at small-scales ~ £, with the assumption 
I <C L p}. Performing this scale separation we obtain, 
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FIG. 1: Time evolution of the large scale magnetic field in dimensionless units in different simulations. The system starting 
from oscillatory behavior, undergoes the transition to reversals regime and finally reaches a steady dynamo for increasing values 
of the Reynolds numbers. 



in the induction equation at large scale, a term which 
describes the action of small scales on the large one con- 
sisting in a turbulent e.m.f. that can be written in terms 
of the Fourier modes of velocity (u(k, t)) and magnetic 
held (b(k, t)) small scale fluctuations as follow: 

e = -^u(k,i) xb*(M) ; (1) 

k 

this is a correlation between velocity and magnetic held 
fluctuations at small scales. Introducing a basis in the 
spectral space: ei(k) , e~2(k) — e3(k)xel(k) , ^(k) = 
ik/|k|; and writing expression ([T]) in a form symmetric 
with respect to the change of k in — k we finally find: 

e = - % [(«* 62 - " 2 & i) + ( M 2 h - ux b* 2 )} (2) 

k(k z >0) 

where u\ and u 2 {b\ and b 2 ) are the components of 
u(k, t) (b(k, i)), along e\ and e 2 . 

We describe the dynamics of the system at large scale 
by integrating the induction equation in local approxima- 
tion: we approximate the toroidal (e^) and the poloidal 
(e~ p ) unit vectors with the cartesian unit vectors, respec- 
tively e x and e z . Hence the field at large scale are: 

u = V(y, z) e x , b = B^{y,t)e x + B p (y,t)e z (3) 

where and B p are the toroidal and poloidal com- 
ponent of the magnetic field, respectively. We describe 
the dynamics at small scales by a shell model [8|, where 
the wave vectors space, in which one considers the MHD 
equations, is divided into a finite number N of shells of 
radius k n = 2 n ko (with n= 0,I,...,N; and fco ~ 2ir/£ is 
the fundamental wave vector). In each shell is assigned 
complex scalar variables u n (t) and b n (t), describing the 
dynamics of velocity and magnetic Fourier modes in the 
shell of wave vectors between k n and fc n+1 . The nonlin- 
ear coupling of neighbor shells is chosen by preserving 



total energy, cross helicity, and magnetic helicity, and, 
at variance with the original MHD shell model, avoiding 
unphysical correlations of phases. We can write the set of 
self-consistent equations for our dynamo model in which 
the e.m.f. is in a form consistent with the shell model: 



dS* * ^ V B, 
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dB p 
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= k n (B^ + B p )b rr 



L 2 ' 



(4a) 
(4b) 
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(4c) 

K+l U n+2) 



—/ik n b n ', 



(4d) 



where the spatial derivative associated at large scale is 
estimated dividing by the typical large scale L; v is the 
viscosity and n is the diffusivity of the MHD flow; n is 
the shell number (n — 0,...N); finally /„ is an exter- 
nal forcing term applied on the first shell (n = 0). This 
forcing, that drives the fluid flow to become unstable, 
has the property to preserve the energy flux to small 
scales. This is the same forcing used in Q with corre- 
lation time equal to 1. The first terms in the RHS of 
Eqs. pc| - pd| describe the effect of the large scale mag- 
netic field on the small-scale turbulent dynamics. Due to 
these terms, the turbulence becomes anisotropic because 
the field fluctuations self organize themselves in coherent 
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FIG. 2: (Color online) Dynamical behavior of Sb/8u for in- 
creasing Rm. Starting from no magnetic field state, we check 
the values of Sb/Su during the running time of the simulation 
in which we increase Rm for v — 10"°. During this time the 
system undergoes dynamo transition (magenta circles) reach- 
ing oscillatory regime, afterwards the system goes towards re- 
versals regime (black triangles) coming to steady states where 
the large scale magnetic field saturates (gray squares). 



structures, channeled along the large scale magnetic field 
bo as Alfven waves. 

We have obtain the e.m.f. consistently with shell model 
without linear approximation. The e.m.f. tends to vanish 
when the system evolves towards a state of strong corre- 
lation between velocity and magnetic field, which char- 
acterizes the Alfvenic subspaces. This is an attractor of 
the dynamics of the system [lOj . We solve the model Eqs. 
assuming V = 0, that is equivalent to solve the dynamo 
problem for Rossby number Ro = 3> 1- Hence 

B p (t) = B^it) = B(t) and the model describes a 2 dy- 
namo problem, in which the shear due to the differential 
rotation at large scale is neglected. Even in the absence 
of the macroscopic shear, the a effect can give rise to 
dynamo action. This effect seems to play a decisive role 
for planetary magnetic fields as geomagnetic field. It can 
be very strong in the fully convective stars as the late- 
type stars and it may provide a possible mechanism to 
explain magnetic activity, along with nonaxisymmetric 
field as observed in many active stars [ill ]. 

The model Eqs. are numerically solved by using a 
fourth order Runge-Kutta scheme. The results are in di- 
mensionless units: the field fluctuations are measured in 
Alfven velocity unit ca, the time in eddy-turn-over time 
(l/(fcoMo))) the lengths are normalized to 1/fco, and fi- 
nally the dissipative coefficients are normalized to c^/fc . 

At the beginning of each simulation, we let the system 
become turbulent at small scales, i.e. we keep B — up 
to 1000 unit times. During this time the kinetic and the 
magnetic energy at small scales grow forming a power law 
spectrum. The system tends to the energy equipartition 
at small scales, exactly achieved when v = [i. After that, 
when we are sure that the turbulence is fully developed, 
we introduce a magnetic field seed of amplitude 10" 10 
at large scale and we check whether the dynamo effect 
starts to develop. 

The numerical results reveal a strong sensitivity of the 
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FIG. 3: The stability curve Rm c vs Pm" 1 in log scale. The 
inset in semi-log scale shows the slight increase of Rm c for 
increasing Pm -1 > 1. 



system with respect to the magnetic Reynolds number 
Rm and a dependence on the hydrodynamic Reynolds 
number Re. Depending on these parameters, the sys- 
tem evolves towards different scenarios: i) no dynamo; 
ii) oscillatory dynamo; iii) magnetic reversals; iv) steady 
dynamo. The transition from oscillatory dynamo to a 
steady dynamo, going through a reversals regime, is ob- 
tained by increasing Re and Rm (Fig.[T]), reproducing in 
this way the same dynamo regimes as observed in the 
liquid metal experiments [l|. 

In order to better characterize the different regimes, 
we have performed a simulation with a fixed value of the 
kinematic viscosity but changing the value of the mag- 
netic diffusivity after a time interval sufficiently long with 
respect to the dynamical evolution times and we have 
studied how the two parameters Sb/du and rms(B)/5u 
change with the magnetic Reynolds number Rm. 

The Su/Sb function illustrated in Fig. [5] for increasing 
Rm shows a zero order discontinuity corresponding to 
dynamo transition (the step for Rm = 67) and a first 
order discontinuity corresponding to the transition from 
oscillatory dynamo to reversals regime. The kind of these 
discontinuities marks the different nature of the transi- 
tions, in particular we can argue that the dynamo effect 
occurs as a plasma instability while the others transitions 
occur with continuity. 

Our model allows us to study the dynamo instability 
inside a very large values of critical parameters not ac- 
cessible by DNS yet, but more realistic for the astrophys- 
ical objects. We have investigated the dynamo threshold 
for different Prandlt number, reproducing the stability 
curve in a very wide range of values, illustrated in Fig. [31 
which shows a little dependence of Rm c for Pm" 1 > 1; 
and a strong monotonic increase of Rm c for decreasing 
Pm" 1 < 1. The reason of this strong increase could be 
related to the need to have small magnetic diffusivity 
to compensate the large viscosity in order to have small 
scales field fluctuations large enough, to consistently trig- 
ger the e.m.f. at large scale for the dynamo onset. On one 
hand, this result is consistent with what was obtained by 
other models in a smaller range of critical parameters: 
Pm" 1 G [1, 100]) [lH, giving us the possibility to be con- 
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FIG. 4: (Color online) Hysteresis cycles of rms(B)/<5« (top 
panel) and Sb/Su (bottom panel) obtained around the stabil- 
ity margin of the bifurcation for the simulation with v = 10~ 
and changing Rm: we increase Rm reaching dynamo onset 
in which the magnetic field displays oscillations (first phase: 
magenta circles), afterwards we decrease Rm (second phase: 
black triangles). 



fident on our model, and on the other hand, it extends 
these results in larger range of values not investigated 
yet. The model hence reaches values more realistic for 
the astrophysical objects. For instance, for interstellar 
or intracluster medium Pm 3> 1, in the Sun's convective 
zone, Pm - 1CT 7 to - 1CT 4 , in planets Pm - 1(T 5 and 
in protostellar disks Pm <c 1 • 

The dynamo onset could take place in subcritical or 
supercritical way. In some models and experiments it 
is subcritical, as it has showed in some small-scale dy- 
namos, revealing a hysteretic behavior fl3l j . In order to 
investigate the nature of this instability as described in 
our model, we have reproduced hysteresis cycles realized 
as following: for a fixed value of is, starting from a state 
of quasi zero-magnetic field characterized by low Rm, 
we have dynamically increased Rm during the simulation 
time in order to destabilize the system from the initial 
state. During the simulation we have checked the time 
evolution of the order parameters: the ratio between the 
r.m.s. of B over Su and 5b over Su. Only after a critical 
value of Rm, corresponding to the threshold for the dy- 
namo action to set up, we have found the increase of the 
large scale magnetic field (first phase) . Once the dynamo 
sets in, we have decreased Rm (second phase). In the sec- 
ond phase, order parameters not follow the path inverse 
to the first phase, realizing a hysteresis loop. We have 
investigated hysteresis cycles considering a wide range of 
values of v down to 10~ n . For the example reported in 
Fig. 2] where v = 10~ 5 , starting from a quasi null mag- 
netic field, we have increased Rm up to Rm c = 67, cor- 



responding to the dynamo threshold, by finding for the 
order parameter the values shown by the magenta cir- 
cles in Fig. |U After the dynamo onset, keep increasing 
Rm, the system displays magnetic oscillation and subse- 
quently magnetic reversals. Once reached these regimes, 
we have decreased Rm. The values of the order param- 
eters, representing now by the black triangles, lie on the 
same curve when the system is in latter regimes, after- 
wards the order parameters lie on a different curve start- 
ing from Rm c to lower values. This shows the hystereti- 
cal behavior of the transition and its subcritical nature. 
We have found also a slight suppression of hydrodynamic 
turbulence during the transition [l4| , revealed by the de- 
crease of hydrodynamic Reynolds number Re ~ Su/kov. 
In the example in Fig. 2J Re, starting from a value around 
8000, slightly decreases to 7000 when the dynamo action 
is developed. We can argue that the subcritical nature of 
dynamo bifurcation could come out from the reduction 
of turbulence due to the magnetic field growing, which 
ensure the maintain of the dynamo effect for lower Rm 
values. 

We present a dynamo model developed for the large 
scale magnetic field. The model solves the induction 
equation in local approximation at large scale, while the 
turbulent dynamics at small scales is described by a MHD 
shell model. We point out that there is no prescribed 
form (like a-term) for e.m.f., but it is consistently de- 
scribed by the analytical form of the shell model. More- 
over at variance with others dynamo shell models [tl [l5[ 
which use ad hoc terms to mimic the dynamo instability, 
this model is able to reproduce different regimes in dy- 
namical way without any ad hoc assumption. An impor- 
tant result of this research is related to the critical role 
that the turbulence plays in the dynamo phenomenon 
showing how different regimes are obtained depending 
on Reynolds numbers of the plasma flow. This property 
of the model is also confirmed by laboratory experiments 
inside the Reynolds number accessible to these. One of 
the strength of our model, due to the use of the shell 
technique, consists of the capability to describe the dy- 
namo transition and the different regimes where it can 
act in a very large values of the Reynolds numbers. The 
study on the dynamo transition shows its subcricital na- 
ture supported by hysteresis cycles around the stability 
margin of the bifurcation. The model reproduces results 
consistent with other dynamo models in the parameter 
range accessible to these, therefore showing its reliability. 
Moreover it makes a step forward in our understanding 
of the dynamo problem due to its capability to cover a 
larger range of values not yet accessible by others mod- 
els and by DNS, but more realistic for a natural dynamo 
processes. Finally we consider the results here obtained 
a good stimulus for both DNS and experimental studies 
on this problem. 
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